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One of the major predictions of inflation is the existence of a stochastic background of cosmo- 
logical gravitational waves (GW). These gravitational waves can induce significant temperature 
anisotropies in the Cosmic Microwave Background (CMB) on the angular scales recently probed 
by the Archeops experiment. Here, we perform a combined analysis of Archeops together with 
information from other CMB experiments and/or cosmological datasets, in order to constrain 
the amplitude of the GW background. We find that, for a scale-invariant GW background, 
the ratio of tensor/scalar perturbations at the CMB quadrupole is now constrained to be 
r < 0.43 at 95% c.L, while the bound on the spectral index of primordial density fluctua- 
tions is n s = 0.97l°;i°- We discuss the implications for future GW detections through CMB 
polarization measurements. 



I. INTRODUCTION 

The last years have seen spectacular advances in our 
ability to confront the inflationary scenario of struc- 
ture formation to observational data. The "multi- 
ple peaks" observed in the Cosmic Microwave Back- 
ground (CMB) angular power spectrum ( [28], [17], 
[23], [32], [35]) are indeed providing strong sup- 
porting evidence for the inflationary predictions of 
a flat universe and of a primordial background of 
scale- invariant adiabatic perturbations (see e.g. [39], 
[29]). More recently, the new CMB results from the 
Archeops experiment ( [1]) have confirmed and re- 
fined the present observational status, sampling angu- 
lar scales between those probed by the COBE satellite 
and the latest high precision datasets. Again, flat- 
ness, adiabaticity and scale invariance are in agree- 
ment with the data ( [2]). 

It has been argued that the next and probably 
most conclusive evidence for inflation would be the de- 
tection of a stochastic background of Gravity Waves 
(GW) (see e.g. [7], [44]). Two types of spacetime met- 
ric fluctuations are indeed naturally produced during 
inflation: density perturbations {scalar modes), which 
form the "seeds" of structure formation, and gravity 
waves {tensor modes) ( [16]). 

The GW background, if detected, would also pro- 
vide valuable information on the inflationary scenario. 
In particular, in most inflationary models (and cer- 
tainly in the simplest ones), the amplitude of the GW 
background is proportional to the square of the energy 
scale of inflation (see e.g. [8]). Furthermore, a comple- 
mentary measurement of the 'tilt' of the GW pertur- 
bations (and of the scalar as well) can give direct in- 
formation up to the second derivatives of the inflaton 
potential, sheding light on the physics at ~ 10 16 GeV 
(see e.g. [19]). 

The GW background leaves an imprint on the CMB 
anisotropies at large scales through the Sachs- Wolfe 
effect. On scales smaller than the horizon at recom- 



bination, however, unlike the anisotropies generated 
by scalar fluctuations, those generated by GW damp 
like fluctuations in a fluid of massless bosons (see e.g. 
[11]). Since the theoretical spectrum, normalized to 
COBE, is a linear sum of the scalar and tensor compo- 
nents, if there is a relevant contribution from GW this 
would lower the predicted amplitude of the acoustic 
peaks on sub-degree angular scales. 

With the advent of the new CMB peaks detections, 
many authors have therefore addressed the question 
of the GW's contribution (see e.g. [26], [21], [39], [12], 
[24], [40], [18]). However, despite the different scale 
dependence, robust constraints on tensor modes re- 
main difficult to obtain. The decrease in the ampli- 
tude of the acoustic oscillations induced by GW can 
indeed be compensated by an increase in one of the 
unconstrained parameters of the model, like, for ex- 
ample, the spectral index of scalar fluctuations ns- 
Therefore, some form of 'cosmic degeneracy' arises in 
the tradeoff between these two (and more) parameters 
(see [26], [12]) and only weak constraints on the GW 
background were obtained. 

In this context, and before more accurate polariza- 
tion data become available (see discussion below) , the 
new results on intermediate angular scales, as recently 
provided by Archeops, can offer an interesting oppor- 
tunity. 

As we illustrate in Fig.l, this spectral region has a 
particular sensitivity to a GW contribution. In the 
figure, we plot two theoretical power spectra. The 
models have identical power on sub-degree scales and 
on COBE scales (considering cosmic variance), but 
different tensor contributions, parametrized by a ten- 
sor over scalar ratio of the angular power spectrum 
quadrupole r = Cj/Cf (see e.g. [21]). 

As we can see, while the two models are degenerate 
on scales I > 200, the degeneracy is broken on larger 
angular scales (see the bottom panel), mostly in the 
region sampled by Archeops. Both increasing n<? and 
adding tensors change the rate of growth of the scalar 
modes from the Sachs- Wolfe plateau towards the first 
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FIG. 1. Best-fit models to recent CMB data with and 
without GW contribution (Top Panel). The Archeops data 
points are shown as open circles. In the Bottom panel we 
plot the % difference between the two degenerate mod- 
els together with the cosmic variance limit (dashed line) 
averaged in bins of A£ = 10. 

peak and this can in principle be used to constrain the 
GW background. 

It is therefore extremely timely to analyze the 
Archeops data allowing the possibility of a GW con- 
tribution in order to see if the amplitude of this back- 
ground can now be better constrained than in the 
past. 

Furthermore, the GW background produces a 
unique statistical signature in the polarization of the 
CMB by inducing a curl component ( [33], [20]), of- 
ten defined as B mode, while scalar (but also ten- 
sor) perturbations produces a gradient component {E 
mode). Given the large number of future and ongo- 
ing CMB polarization experiments, it is interesting to 
forecast from the present CMB temperature data the 
expected amplitude of the B modes and/or if the E 
modes produced by tensors can be distinguished from 
those produced by scalar perturbations only. 

We pursue this investigation in the present Rapid 
Communication as follows: in Section II we illustrate 
our analysis method. In section III we present our 
results. Finally, in section IV, we discuss our findings. 

II. ANALYSIS: METHOD 

As a first step, we consider a template of flat, adi- 
abatic, A-CDM scalar and tensor spectra computed 
with CMBFAST ( [36]), sampling the various param- 
eters as follows: il c d m h 2 = u c d m = 0.05, ...0.25, in 
steps of 0.02; n b h 2 = u b = 0.009, 0.024, in steps of 



0.003, fl A = 0.5, ...,0.95, in steps of 0.05. Our choice 
of the above parameters is motivated by the Big Bang 
Nucleosynthesis bounds on u>t, (both from D [6] and 
A He + 7 Li [9]), from supernovae ( [14]) and galaxy 
clustering observations (see e.g. [38]). 

Variations in the tensor and scalar spectral indices, 
ns and n-r are not computationally relevant. How- 
ever, we restrict our analysis to relevant inflationary 
values ns — 0.7, 1.3 and we fix nx = (see discus- 
sion below for different values of Ut)- 

Furthermore, the value of the Hubble constant 
is not an independent parameter, since h = 
\J {^cdm + w f))/(l — ^a)- We also include the further 
top-hat prior h = 0.7 ± 0.2 ( [13]) and we consider 
only models with age t n > 11 Gyrs. 

We allow for a reionization of the intergalactic 
medium by varying the compton optical depth pa- 
rameter t c in the range r c = 0.0, ...,0.45 in steps of 

0. 05. We note here that high values of r c are in severe 
disagreement with recent estimates of the redshift of 
reionization z re ~ 6 ± 1 (see e.g. [15]) which points 
towards r c ~ 0.05 — 0.10. On the other hand, if the 
reported CBI excess at £ ~ 3000 is due to Sunyaev- 
Zeldovich effect, then this would favour values r c ~ 0.3 
( [3])- 

For the CMB data, we use the recent results from 
the BOOMERanG-98, DASI, MAXIMA-1, CBI, VSA 
and Archeops experiments. The power spectra from 
these experiments were estimated in 19, 9, 13, 14, 10 
and 16 bins respectively (for the CBI, we use the data 
from the MOSAIC configuration, [10]), spanning the 
range 2 < £ < 1500. We also use the COBE data from 
the RADPACK compilation ( [34]). 

For the CBI, DASI, MAXIMA-I and VSA exper- 
iments we use the publicly available correlation ma- 
trices and window functions. For the Archeops and 
BOOMERanG experiments we assign a flat interpola- 
tion for the spectrum in each bin £(£ + l)Ce/2-K — Cb, 
and we approximate the signal Cb inside the bin to 
be a Gaussian variable. The likelihood for a given 
theoretical model is defined by — 21n£ = (C 1 ^ — 
C B X )M B B'{C%, - Cf,) where M BB < is the Gaussian 
curvature of the likelihood matrix at the peak. 

We consider 5%, 10%, 4%, 5%, 3.5% and 5% 
Gaussian distributed calibration errors (in AT) for 
the Archeops, BOOMERanG-98, DASI, MAXIMA- 

1, VSA, and CBI experiments respectively and we 
include the beam uncertainties by the analytical 
marginalization method presented in ( [4]). 

Finally, we parametrize the GW contribution by the 
tensor over scalar quadrupole ratio r = Cj/Cf and 
we rescale the sum spectrum by a prefactor Cio, as- 
sumed to be a free parameter, in units of C^f BE . 

III. ANALYSIS: RESULTS 

The main results of our analysis are plotted in Fig. 2. 
In the left top panel we plot the likelihood contours in 
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FIG. 2. 68%, 95% and 99% confidence regions in the 
ns —r (Top Panel, Left), ns — t (Top Panel, Right), r — r 
(Bottom Panel) planes for the models considered in our 
analysis (see text). The line contours are confidence levels 
without the Archeops data. 



the ns — r plane, maximizing over the remaining nui- 
sance parameters. As we can see, in the framework of 
models we considered, the gravitational wave contri- 
bution is constrained to be r < 0.2 (r < 0.43) at 68% 
C.L. (95% C.L.), with n s = 0.97™ (68% C.L.). 
While the inclusion of the Archeops data has little ef- 
fect on ns, it drastically improves the constraint on 
r. Removing the Archeops data yields r < 0.6 at 95% 
C.L.. 

In the right top panel of Fig. 2, we plot the likeli- 
hood contours in the ns — t c plane. As we can see, 
the present CMB constraint on r c is rather weak, with 
r c < 0.25 (r c < 0.36) at 68% C.L. (95% C.L.). It is in- 
teresting to note that the inclusion of the ARCHEOPS 
datapoints has little effect. 

Finally, in the bottom panel of Fig. 2, we plot 
the likelihood contours in the r — t c plane. An in- 
crease in r c or r produces a similar damping on the 
small/intermediate angular scales. It is interesting to 
notice that the present data is allowing just a well 
defined amount of small-scale damping. Values of 
r c ~ 0.3 are in disagreement with the presence of a 
tensor component. If r c > 0.2 then r < 0.05 at 68% 
C.L.. 

To each theoretical model in the likelihood planes 
produced in Fig. 2, is possible to associate a theoretical 
polarization power spectrum and translate the confi- 
dence contours into an expected maximum and min- 
ima polarization signal. 

We do this in the 3 panels of Fig. 3, where we plot 
the envelope of the minima and maxima polarization 
spectra that, in the panels of Fig. 2, are at 95% c.l. 
consistent with the CMB temperature data. 
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FIG. 3. Maximum and minimum levels of tempera- 
ture-polarization cross correlation (Top Panel), B-modes 
(Central Panel) , E- modes (Bottom Panel) allowed at 95% 
C.L. from present CMB temperature data under the as- 
sumption of the models described in the text. 



As we can see from the center panel of Figure 3, 
the level of the -B-modes, is expected to be of ~ 0.2 
/iK, at maximum. The signal is out of the reach of 
most of the current polarization experiment like DASI 
or POLATRON which are sensitive to few fiK. Near 
future experiments like B2K or QUEST, will probably 
have enough sensitivity to have a statistical B-mode 
detection. However, the S-signal in the angular re- 
gion sampled by these experiments (£ > 50), can be 
contaminated by a foreground component due to the 
conversion of E modes to B modes from gravitational 
lensing (see Fig. 3) ( [41]). Higher-order correlations 
will be necessary to map the cosmic shear and sub- 
tract this contribution to the B mode ( [30]). 

Tensor perturbations produce E modes as well. 
However, the amplitude of the E tensor modes is pre- 
dicted to be generally much smaller than those from 
the scalar modes (see bottom panel) . A window of op- 
portunity may appear in the temperature-polarization 
(< TE >) cross-correlation spectra, where, at I ~ 50, 
the amplitude from tensor can be larger than those 
from scalar modes, leaving a possible detectable ex- 
cess for experiments like QUEST or B2K. 

In order to cross-check if any information can be 
obtained on riT we performed the analysis on just one 
cosmological model defined by £1a = 0.7, uJb = 0.022, 
Qtot = 1, t c = 0.04. We then considered ten- 
sor contributions by varying the scalar and tensor 
spectral indices independently: ns = 0.7, ...,1.3 and 
n T = -0.3, ...,0.0, step 0.01. We found that the ten- 
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sor spectral index is not constrained by the present 
data, but that a value of tit = is preferred. 



IV. CONCLUSIONS 

In this Rapid Communication we have presented 
new constraints on the stochastic background of gravi- 
tational waves from recent microwave anisotropy data. 
Thanks to Archeops, our results improve the con- 
straints on tensor modes from previous analyses (see 
e.g. [40], [21]). 

In the framework of models we considered, we found 
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(at 95% C.L.) r < 0.43 and n s = 0.97 
energy scale of inflation E\ n f can be related to tensor 
by Ef n j = 0.65Cjmp Z . The above bound translates 
therefore in E inf < 1.6 x W 16 GeV. 

When comparing with the results presented in [2] , a 
part from the different template of theoretical models 
considered, our analysis differs mainly in the following 
points: we assumed the low-£ Archeops bins as gaus- 
sian distributed, we included the COBE data using 
the RADPACK compilation, we have a strong upper 
limit on lui, < 0.025 from BBN and, finally, we numer- 
ically computed the models with t c > (while in [2] 
an analytical formula was used). 

The GW background induces a unique signature in 
the polarization of the CMB by producing a curl com- 
ponent, not present in the case of scalar perturbations. 
In the set of models we considered (and under the as- 
sumption of a bayesian method of statistical analy- 
sis) we found that the maximum expected level of B 
modes allowed by current data is of about <~ 0.2[iK, 
which can be partially attainable by near future exper- 
iments and severly contaminated by lensing E — > B 
conversion. 

The E modes expected from gravity waves are lower 
than the E modes expected from scalar perturbations. 
However, the tensor < TE > cross-correlation might 
be larger at I ~ 50. 

In this context, deviations in the < TE > cross- 
correlation scalar spectrum at I ~ 50 can possibly 
offer competitive information with respect to B modes 
search. 

All our predictions derived from the temperature 
data are consistent with the recent claim of detection 
of polarization E modes from the DASI experiment ( 
[22]). 

As a final remark, we want to stress that the results 
presented here have been obtained under the assump- 
tion of a theoretical framework. The bounds on the 
polarization spectrum must be considered just as an 
indication of what future observations may detect. 

In particular, we just considered ut = in the 
main analysis and we looked at the effect of having 
ut as low as —0.3. For inflation, only values Ut < 
can be considered and we checked that varying ut 
has little effect on the final results. Spectra with 
"blue" (nr > 0) spectral indices can be produced in 



Ekpyrotic ( [37]) or Prc-Big Bang (see e.g. [27]) sce- 
narios. However, extremely blue spectra (nr ~ 2) 
are excluded by constraints on the GW energy den- 
sity background from timing milli-second binary pul- 
sars [31]. Allowing for extra primordial perturbation 
modes like isocurvature, will probably tight our con- 
straints on GW, since the shape of CDM scalar isocur- 
vature modes is similar to those from adiabatic tensor 
modes. However, considering the most general initial 
conditions scheme and including cross correlations, 
will certainly enlarge our constraints ( [5]). Includ- 
ing curvature (Cltot ^ 1) would relax our bounds on r 
(see e.g. [39]). Non-flat models in agreement with the 
CMB data arc in general closed models, which, a part 
from a few exceptions ( [25]), arc difficult to obtain 
from inflation. Finally, including a different model for 
dark energy like quintessence would change the large 
scale anisotropy through the Integrated Sach Wolfe 
effect (sec e.g. [47]), affecting our constraints as well. 

Even if the results presented here do not hint for a 
presence of GW background, the data is still consis- 
tent with a sizable tensor contribution. It will there- 
fore be the duty of future and ongoing experiments to 
scrutinize this fundamental prediction from inflation. 
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